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Abstract: The upcrossings index 0 < ry < 1, a measure of the degree of local dependence in the up¬ 
crossings of a high level by a stationary process, plays, together with the extremal index 9, an important 
role in extreme events modelling. For stationary processes, verifying a long range dependence condition, 
upcrossings of high thresholds in different blocks can be assumed asymptotically independent and there¬ 
fore blocks estimators for the upcrossings index can be easily constructed using disjoint blocks. In this 
paper we focus on the estimation of the upcrossings index via the blocks method and properties such 
as consistency and asymptotic normality are studied. We also enlarge the family of runs estimators of 
r] and provide an empirical way of checking local dependence conditions that control the clustering of 
upcrossings to improve the estimates obtained with the runs method. 

We compare the performance of a range of different estimators for rj and illustrate the methods using 
simulated data and financial data. 
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1 Introduction 

Extreme Value Theory aims to predict occurrence of rare events but with disastrous impact, making 
an adequate estimation of the parameters related with such events of primordial importance in Statistics 
of Extremes. The extremal index 9 G (0,1] and the upcrossings index rj G (0,1] play an important role 
when modelling extreme events. The knowledge of these parameters, 9 and ij, entails in particular the 
understanding of the way in which exceedances and upcrossings of high levels, respectively, cluster in 
time. They provide different but complementary information concerning the grouping characteristics of 
rare events. 

Suppose that we have n observations from a strictly stationary process X = {Xn}n>i, with marginal 
distribution function F and finite or infinite right endpoint xp = supjx ; F{x) < 1}. X is said to have 
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( 1 . 1 ) 


extremal index 9 G [0,1] if for each r > 0 there exists a sequence of thresholds {un^}n>i such that 


n(l — F{u^^^)) ->• r and Fn{u^n^) = P{Mn < ->• e 

n—)-H-oo n^+oo 

where = max{Xi,..., (Leadbetter (1983) [E]). This parameter is a measure of clustering 
tendency of extremes, more precisely 6~^ is the limiting mean cluster size in the point process of ex¬ 
ceedance times over a high threshold, under an appropriate long range dependence condition (Hsing et 
al. (1988) [H]). For independent and identically distributed (i.i.d.) processes, extremes are isolated and 
hence 0 = 1, whereas for stationary processes, the stronger the dependence the larger the clusters of 
exceedances and therefore the smaller 9. 

Many results in extreme value theory may be naturally discussed in terms of point processes. When 
looking at the point process of upcrossing times of a threshold Un, that is 


Nn{Un){B) — ^{Xi<Un<Xi+i}^i. {B), B C (0,1], 


( 1 . 2 ) 


2=1 


where Ta denotes the indicator of event A and 5a the unit mass at a, Ferreira (2006) [6] showed that, under 
a long range dependence condition, if the sequence of point processes of upcrossings Nn{un) converges in 
distribution (as a point process on (0,1]), then the limit is necessarily a compound Poisson process and 
the Poisson rate of the limiting point process is rjv, when the limiting mean number of upcrossings of 
is > 0. That is, there is a clustering of upcrossings of high levels, where the underlying Poisson values 
represent cluster positions and the multiplicities are the cluster sizes. In such a situation there appears 
the upcrossings index r/ which has an important role in obtaining the sizes of the clusters of upcrossings. 

The stationary process X is said to have upcrossings index r/ G [0,1] if for each u > 0 there exists a 
sequence of thresholds {u^^}n>i such that 


nP{Xi < ul^'> < X 2 ) -> and P(X„(m^^))(( 0, Ij) = 0) -^ e-^r 


(1.3) 


71^+00 n^+oo 

The upcrossings index 77 exists if and only if there exists the extremal index 9 (Ferreira (2006) [ 6 ]) and 
in this case 

T] = - 9. (1.4) 

Note that for a high level Un such that nP{Xi > Un) = t' and nP{Xi < Un < X 2 ) = F we have 
r' > 17 ' > 0 and therefore rj >9. 

Let us now consider the moving maxima sequence of Ferreira (2006) [6] in the following example: 

Example 1.1 Let {Yn}n >-2 be a sequence of independent uniformly distributed on [0,1] variables with 
common distribution function F. Dehne the 3-dependent moving maxima sequence as 
Xn = max{Yn,Yn- 2 ,yn- 3 }, n > 1, and consider {Zn}n>i a sequence of i.i.d. random variables with 
common distribution F^. The underlying distribution of X^ is also F^ and it has an extremal index 


0 = 1/3 and upcrossings index 77 = 1/2. Moreover, for «„ = 1 — r > 0, it holds nP{Xi > Un) 
and nP{Xi < Un < X 2 ) 


n—>-+oo 


-G 3r 


n—>-+oo 


> 2r, agreeing with (1.4). In Figure 1 we can see the sizes equal to 3 and 


2, respectively for the clusters of exceedances and upcrossings of a high level by the max-autoregressive 
sequence Xn- Note that Xn verihes condition D^^\un) of Chernick et al. (1991) [3] (see Ferreira (2006) 
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El), therefore two runs of exceedances separated by one single non-exceedance are considered in the same 
cluster. 

A shrinkage of the largest and smallest observations for the 3-dependent sequence can in general also 
be observed, despite the fact that we have the same model underlying both sequences. 



Figure 1: Sample paths from an i.i.d. process and a 3-dependent process. The horizontal line corresponds 
the the 80th percentile of n = 50 realizations of the 3-dependent sequence n > 1. 


In the i.i.d. setting we obviously have an upcrossings index r/ = 1. This also holds for processes 
satisfying the local dependence condition D"{un) of Leadbetter and Nandagopalan (1989) [T7|, which 
enables clnstering of upcrossings. Other properties of the upcrossings index can be found in Ferreira 
(2006, 2007) [6], [7] and Sebastiao et al. (2010) [23], relations with upcrossings-tail dependence coefficients 
are presented in Ferreira and Ferreira (2012) |8|. 

It is clear that clustering of extremes can have catastrophic consequences and upcrossings of a high 
threshold can also lead to adverse situations. These motivate an increasing interest in the grouping 
characteristics of rare events. Applications include, for example, the temporal distribution of large 
financial crashes or the evaluation of oscillation periods in financial markets. A search for reliable tools 
to describe these features has become constant in the past years, a key aspect being the estimation of the 
extremal index 9 and the npcrossings index rj which govern, respectively, the clustering of exceedances 
and upcrossings of a high level by an univariate observational series. 

Several estimators for the extremal index 9 can be found in the literature (see Ancona-Navarrete and 
Tawn (2000) [I], Robert et al. (2009) |2T] and references therein). Most of the proposed estimators 
are constructed by the blocks method or by the runs method. These methods identify clusters and 
construct estimates for 9 based on these clusters. The way the clusters are identified distinguishes both 
methods. The blocks estimators are usually constructed by using disjoint blocks because exceedances over 
high thresholds of different blocks can be assumed asymptotically independent. More recently, an inter¬ 
exceedance times method has been proposed that obviates the need for a cluster identification scheme 
parameter (Ferro and Segers (2003) |9|). In what concerns the estimation of the upcrossings index, little 
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has yet been done. In Sebastiao et al. (2013) [Mj we find a runs estimator and a blocks estimator 
derived from a disjoint blocks estimator for the extremal index, given in m, for this parameter. It was 
shown that the runs estimator has a smaller asymptotic variance and a better performance, nevertheless 
it assumes the validation of a local dependence condition which does not hold for many well known 
processes and can be cumbersome to verify. 

This paper focuses on the estimation of r]. The main novelty is our proposal of new blocks estimators 
and the study of their properties, namely consistency and asymptotic normality, as well as an improvement 
of the existing runs estimator. We examine the behaviour of the new and existing estimators of t] and 
assess their performance for a range of different processes. Assuming the validation of local dependence 
conditions D^^\un), k > 3, given in Ferreira (2006) [B] we define new runs estimators for rj. With these 
new estimators we obtain better estimates of r], comparatively to the ones obtained in Sebastiao et al. 
(2013) [23] where k = 3 was considered. We analyze a way of checking conditions D^^\un) preliminary 
to estimation, which helps considerably when using the runs estimators to estimate rj. 

We consider the problem of estimating the upcrossings index in Section 2. There we suggest an 
estimator motivated by the relation 


+O0 


E[w„(nM) I > 0] = 


n—>-+00 


> r] 


-1 


(1.5) 


where 


:— P(yN.p^{Un) — j I k^rni.'^n) ^0) 


is the conditional cluster size distribution and Nr^{un) = A'„(u„)((0,r„/n]), under suitable choices of 

rn - ——> + 00 . This relation is a consequence of Lemma 2.1 of Ferreira (2006) [6], that states that 

under condition A(m„) of Hisng et al. (1988) [ 13 ] upcrossings over disjoint blocks are asymptotically 


independent, and definition (1.3) of the upcrossings index rj. It is extremely important in the interpretation 
of the upcrossings index, since it tells us that the reciprocal of the upcrossings index can be viewed as 
the limiting mean cluster size of upcrossings. 

Other blocks estimators are also presented and their properties studied in Section 2. In subsection 
2.2 we revisit the runs declustering method considered in Sebastiao et al. (2013) [ 23 ] and define new runs 
estimators for r] under the validation of any dependence condition D^^\un), k > 3. A way of checking 
conditions D^^\un), k > 3, preliminary to estimation is also given. 

We compare, in Section 3, the performance of the estimators through a simulation study for a range 
of different processes with r] < 1 to rj = 1 and, in Section 4, by the application of the methods to the 
DAX daily log returns. Conclusions are drawn in Section 5 and proofs are given in the appendix. 


2 Estimation of the upcrossings index 

We shall now consider the problem of estimating the upcrossings index in the following setting. As 
previously mentioned X = {Xn}n>i is a stationary sequence of random variables with a continuous 
marginal distribution function F and a nonzero upcrossings index rj. Assume also that condition A{un'^) 
of Hsing et al. (1988) [13] holds for X. 
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2.1 Blocks declustering 

The blocks declustering scheme consists in choosing a block length rn and partitioning the n obser¬ 
vations (Xi,... ,Xn) of X into kn = [n/r„] disjoint blocks, where [y] denotes the integer part of y. Each 
block that contains at least one upcrossing is treated as one cluster. If in Figure 1 we consider blocks of 
size 10 (r„ = 10) we can identify 3 clusters. 


2.1.1 A disjoint blocks estimator 


Under relation (1.5), rj ^ is the (asymptotic) mean cluster size of upcrossings, and so a natural way 


to estimate y is as the reciprocal of the sample average cluster size, that is, count the total number of 
blocks with at least one upcrossing and divide it by the total number of upcrossings of a high threshold 
u. We then propose the following blocks estimator for rj : 

-1 


y^iur 


I +00 




( 2 . 6 ) 


where 




7rn(j; Un) ■— ^ ^ 

i=l / i=l 

is the empirical conditional distribution of Nr'^}{un) and Nr’^}{un) = Nn{un){[{{i — l)rn + 

YlJ=(i-i)r - 1-1 ^{Xj<u„<Xj+i}^ 1 < * < kn, represents the number of upcrossings of the level Un in the 
ith—block of length r^. 


(2.7) 

l)/n, irn/ri]) = 


As we can see from (2.6) and (2.7) the average number of upcrossings in the blocks estimates rj and 


the proportion of blocks with j upcrossings estimates TT{j]Un) = lim„^+oo Xn{j; Ur. 


The level Un in (2.6) is a tuning constant that determines the quality of the estimate. If Un is too 


high there will be just a few upcrossings of the threshold so estimates will be highly variable, but if Un is 
too low the approximation to the limiting characterization may be poor so the estimates may be biased. 
The sample size is important here as it influences what is a too high threshold. 

Note that when considering n finite the sequences {rn}n>i and {kn}n>i determine the cluster iden¬ 
tification. The values of are determined by the dependence structure of the process. If a process has 
weak dependence, the high level upcrossings will occur in well defined isolated groups, so small values 
for r„ will suffice, similarly large values of are required with strong dependence. 

As we shall see in the next section, the quality of the estimate depends heavily on the choice of both 
parameters, Un and r„. 


2.1.1.1 Consistency and asymptotic normality 


Consistency of the estimator can only be achieved with lower thresholds than the ones considered 


in (1.3), since for these levels there are insufficient upcrossings to give statistical “consistency” for the 
estimator. That is, as n increases the value of rjn does not necessarily converge appropriately to the 
value rj. We shall therefore consider thresholds Vn = '^[^/c ]’ some fixed u > Q, that satisfy 


nP{Xi < Vn < X 2 ) ~ Cni^ as n —)■ -|-oo. 


( 2 . 8 ) 
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where {cn}n>i is a sequence of real numbers such that 
0, for a sequence of integers {kn}n>i satisfying kr, 


n^+oo 


n^+oo 
-)■ +00 


> + 00 , n/Cn 


n^+oo 


> +00 and Cn/k„ 


n^+oo 


The following result shows that for thresholds Vn, satisfying (2.8), convergence (1.5) still holds. 


Theorem 2.1 Let {r„}n>i) {kn}n>i and {cn}n>i be sequences of real numbers such that 

Cn —)■ +C) 0 , kn +C) 0 , > + 00 , Cn/kn —> 0, < n knrn ~ n, as +00 


and {u,i}„,>i a sequence of thresholds defined by (2.8). Suppose there exists In = o{rn) such that 


^n.l 


- > 0 , where oini„, n> 1, are the mixing coefficients of the A(u„) condition. 

’ n^+oo ’ 


Then 


lim —P{NrSvn) > 0) = 
n^+00 Cn 


(2.9) 


and (1.5) holds for Vn, that is 


^[^r„iVn) I dSfr„{Vn) > 0] 


n^+cxo 


> 7? 


Note that from Theorem 2.1 if (2.8) holds then the order of magnitude of the expected number of blocks 
having upcrossings of Vn is Cn - > +oo. This justifies the need to consider lower levels satisfying (2.8) 

n^+oo 


in order to guarantee the consistency of the estimator . Such levels were also considered in Sebastiao 
et al. (2013) [23| to obtain the consistency of the runs estimator of ij. 

Following the style of proofs used by Hsing (1991) [15], we show in the next results that the blocks 
estimator of the upcrossings index is a consistent estimator for rj and that it is asymptotically normal. 
These results require assumptions on the limiting behaviour of the first and second moments of Nr„{vn). 


Theorem 2.2 Suppose that the conditions of Theorem 2.1 hold and for some u > 0, 


( 1 ) 


>’n n^+oo 


> 0 , 


( 2 ) 

Then 


n^+oo 


0 . 


hn (Vn) 


n^+oo 


> T]. 


The underlying idea to obtain the asymptotic distribution of the estimator is to split each block 
into a small block of length In and a big block of length rn — ln- Since = o(r„) and u -^ 0 it is 

’ n—>+ 2 x 3 

ensured that In is sufficiently large such that blocks that are not adjacent are asymptotically independent, 
but does not grow to fast such that the contributions of the small blocks are negligible. 


Theorem 2.3 If conditions of Theorem \2.1\ hold and for some u > 0, 

( 1 ) 0, for all e > 0, 

(2) ^E[Nl{vn)] 


n^+oo 


(D„)>ec„}J 

> rjua'^, for some cr^ < + 00 . 
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Then 


5^(iV»(u„)-i?[iV«(u„)]) 


i=l 


- 2=1 




( 2 . 10 ) 


Remark 1 From assumption (2) of the previous theorem we have that E[Nf (vn) \ Nf (vn) > 0] - > 

" " n^+oo 

CT^, for some al^ > 0, which means that Var[N^^{vn) \ N‘^^{vn) > 0] - -—> — ri~‘^, the asymptotic 

variance of a cluster size. 

The asymptotic normality of rj^ is now an immediate consequence of the previous result. 

Corollary 2.1 If the conditions of Theorem 2.3 hold then 

- Vn) —^ AA [o, - 1 )) , 

n—>-+oo \ V J 

where rjn = 

From the previous result it is clear that the asymptotic distribution of is not well defined when 
cr^ < since in this case a smaller or equal to zero variance is obtained. 


Remark 2 The asymptotic distribution found in Corollary 2.1 differs from the one obtained with the 
runs estimator in Sebastido et al. (2013) only on the value of a'^. Thus, the value of determines 
which of the two is the most efficient estimator of rj. 


In order to control the bias of the upcrossings index estimator 7]((, lets assume that the block sizes 
are sufficiently large so that as n —)■ +oo 

V^iVn-v) -^ 0 . ( 2 . 11 ) 

n^+oo 

The asymptotic variance of will depend on p, v and which can easily be estimated with consistent 
estimators as shown in the next result. This result enables the construction of approximate confidence 
intervals or a hypothesis test regarding rj. From a practical viewpoint it is more useful than Corollary 2.1 
and follows readily from the previous results. 


Corollary 2.2 If the conditions of Corollary 2.1 hold, (2.11) and 


k 


r2 „^+oo^ °- 


Then 


I EfeiJvgK) 

VniiVn^n)^ “ 1 ) 


iVn - V) 


n^+oo 


>AA(0,1), 


where 0}{vn) = E7=(i-i)r-„+i ^{Xi<v„<x,+,} and a'i = 


2 _ 




( 2 . 12 ) 
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Given a random sample (Xi,..., Xn) of X and a threshold u we obtain, from Corollary 2.2, the 
asymptotic 100(1 — a) percent two-sided confidence intervals for r] 




Vn ^1—a/2-| 


I VniiVn^n)^ - 1) -S , ^ 


I VniiVn^n)^ - 1 ) 

EfciiviJw , 


(2.13) 


where 2;i_a/2 denotes the 1—a(/2 standard normal quantile, = 


B _ {ivff(»)>0} , ^2 ^ _ 

Z^*=l {jv('bu)>0} 


.2 _ E-^1 


The intervals (2.13) are approximations of the true confidence intervals for finite samples. 


The finite sample properties of are now analyzed in simulated samples from three different processes 
with r] < 1 to r] = 1. The results for all examples were obtained from 5000 simulated sequences of size 
n = 5000, of each process. In each example, the tables contain Monte-Carlo approximations to biases and 
root mean squared errors (RMSE), in brackets, of the estimator (u) and the 95% confidence intervals 


based on the asymptotic Normal approximation for the blocks estimates (2.13), for different values of 
rn and different thresholds corresponding to the range between the 0.7 and 0.99 empirical quantiles, q,. 
Note that the confidence intervals (2.13) can only be computed when > 1/(77,^)^. 


Example 2.1 (MM process) Let Xn = maxjl^, Yn- 2 , Yn-s}, n > 1, denote the moving maxima 
process, where {Yn}n >-2 is a sequence of independent and uniformly distributed on [0,1] variables, 
considered in Example 1.1. 


rn 

90.7 

90.8 

90.9 

90.95 

90.99 

5 

0.1983 ( 0 . 1986 ) 
( 0 . 6775 , 0 . 7191 ) 

0.1956 ( 0 . 1960 ) 
( 0 . 6716 , 0 . 7196 ) 

0.1964 ( 0 . 1972 ) 
( 0 . 6641 , 0 . 7287 ) 

0.1983 ( 0 . 2000 ) 
( 0 . 6536 , 0 . 7430 ) 

0.2001 ( 0 . 2093 ) 
( 0 . 6029 , 0 . 7974 ) 

7 

0.0991 ( 0 . 0997 ) 
( 0 . 5788 , 0 . 6195 ) 

0.1081 ( 0 . 1088 ) 
( 0 . 5850 , 0 . 6312 ) 

0.1228 ( 0 . 1240 ) 
( 0 . 5927 , 0 . 6530 ) 

0.1322 ( 0 . 1345 ) 
( 0 . 5916 , 0 . 6729 ) 

0.1402 ( 0 . 1514 ) 
( 0 . 5546 , 0 . 7258 ) 

10 

- 0.0034 ( 0 . 0102 )+ 
( 0 . 4765 , 0 . 5167 ) 

0.0215 ( 0 . 0244 ) 
( 0 . 4984 , 0 . 5446 ) 

0.0561 ( 0 . 0583 ) 
( 0 . 5266 , 0 . 5856 ) 

0.0770 ( 0 . 0801 ) 
( 0 . 5389 , 0 . 6151 ) 

0.0948 ( 0 . 1069 ) 
( 0 . 5208 , 0 . 6687 ) 

15 

- 0.1154 ( 0 . 1157 ) 
( 0 . 3658 , 0 . 4035 ) 

- 0.0726 ( 0 . 0733 ) 
( 0 . 4042 , 0 . 4506 ) 

- 0.0126 ( 0 . 0196 )+ 
( 0 . 4570 , 0 . 5178 ) 

0.0250 ( 0 . 0326 ) 
( 0 . 4873 , 0 . 5627 ) 

0.0583 ( 0 . 0727 ) 
( 0 . 4960 , 0 . 6207 ) 

25 

- 0.2418 ( 0 . 2419 ) 
( 0 . 2430 , 0 . 2733 ) 

- 0.1882 ( 0 . 1883 ) 
( 0 . 2906 , 0 . 3331 ) 

- 0.0995 ( 0 . 1004 ) 
( 0 . 3688 , 0 . 4322 ) 

- 0.0369 ( 0 . 0420 ) 
( 0 . 4229 , 0 . 5032 ) 

0.0243 ( 0 . 0451 ) 
( 0 . 4699 , 0 . 5787 ) 

50 

- 0.3652 ( 0 . 3652 ) 
( 0 . 1261 , 0 . 1436 ) 

- 0.3228 ( 0 . 3228 ) 
( 0 . 1625 , 0 . 1919 ) 

- 0.2246 ( 0 . 2249 ) 
( 0 . 2461 , 0 . 3047 ) 

- 0.1300 ( 0 . 1314 ) 
( 0 . 3261 , 0 . 4138 ) 

- 0.0136 ( 0 . 0385 ) 
( 0 . 4296 , 0 . 5432 ) 

100 

- 0.4324 ( 0 . 4324 ) 
( 0 . 0632 , 0 . 0719 ) 

- 0.4095 ( 0 . 4095 ) 
( 0 . 0826 , 0 . 0984 ) 

- 0.3403 ( 0 . 3403 ) 
( 0 . 1392 , 0 . 1801 ) 

- 0.2400 ( 0 . 2404 ) 
( 0 . 2193 , 0 . 3007 ) 

- 0.0568 ( 0 . 0693 ) 
( 0 . 3713 , 0 . 5152 ) 

1000 

- 0.4932 ( 0 . 4932 ) 
( 0 . 0064 , 0 . 0071 ) 

- 0.4909 ( 0 . 4909 ) 
( 0 . 0084 , 0 . 0097 ) 

- 0.4836 ( 0 . 4836 ) 
( 0 . 0145 , 0 . 0183 ) 

- 0.4685 ( 0 . 4685 ) 
( 0 . 0260 , 0 . 0369 ) 

- 0.3476 ( 0 . 3478 ) 
( 0 . 0946 , 0 . 2102 ) 


Table 1: Bias and root mean squared error (RMSE), in brackets, of the estimator //^^(u) (first line) and 
the 95% confidence intervals based on the asymptotic Normal approximation for the blocks estimates 


(2.13) (second line), for an MM process with 9 = 1/3 and rj = 1/2. Best values are ticked with a plus. 


Example 2.2 (AR(1) process) Let Xn = —j^Xn-i + (-n: n > 1, denote a negatively correlated au¬ 
toregressive process of order one, where {e„}„>i is a sequence of i.i.d. random variables, such that, for a 
fixed integer /3 > 2, e„ t/{i,...,^,l} andXo [7(0,1) independent of e„. 

Condition D^^\un) holds for this stationary sequence, it has extremal index 9 = 1 — l/fP' and 
upcrossings index — 1/732^ /3 > 2 (see Sebastiao et al. (2010) p3]L 
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Tn 

90.7 

90.8 

90.9 

90.95 

90.99 

5 

- 0.1290 ( 0 . 1291 ) 
( 0 . 6064 , 0 . 6355 ) 

0.0143 ( 0 . 0170 ) 
( 0 . 7438 , 0 . 7848 ) 

0.0881 ( 0 . 0894 ) 
( 0 . 8092 , 0 . 8670 ) 

0.1008 ( 0 . 1030 ) 
( 0 . 8105 , 0 . 8910 ) 

0.1031 ( 0 . 1128 ) 
( 0 . 7659 , 0 . 9402 ) 

7 

- 0.2834 ( 0 . 2834 ) 
( 0 . 4547 , 0 . 4784 ) 

- 0.1245 ( 0 . 1248 ) 
( 0 . 6052 , 0 . 6457 ) 

0.0098 ( 0 . 0190 )+ 
( 0 . 7279 , 0 . 7916 ) 

0.0566 ( 0 . 0611 ) 
( 0 . 7617 , 0 . 8514 ) 

0.0747 ( 0 . 0902 ) 
( 0 . 7281 , 0 . 9213 ) 

10 

- 0.4176 ( 0 . 4176 ) 
( 0 . 3237 , 0 . 3411 ) 

- 0.2732 ( 0 . 2732 ) 
( 0 . 4591 , 0 . 4944 ) 

- 0.0874 ( 0 . 0888 ) 
( 0 . 6290 , 0 . 6962 ) 

0.0009 ( 0 . 0235 )+ 
( 0 . 7022 , 0 . 7996 ) 

0.0497 ( 0 . 0730 ) 
( 0 . 6964 , 0 . 9030 ) 

15 

- 0.5278 ( 0 . 5278 ) 
( 0 . 2165 , 0 . 2279 ) 

- 0.4197 ( 0 . 4197 ) 
( 0 . 3172 , 0 . 3435 ) 

- 0.2139 ( 0 . 2143 ) 
( 0 . 5037 , 0 . 5685 ) 

- 0.0776 ( 0 . 0810 ) 
( 0 . 6211 , 0 . 7237 ) 

0.0220 ( 0 . 0605 ) 
( 0 . 6635 , 0 . 8805 ) 

25 

- 0.6166 ( 0 . 6166 ) 
( 0 . 1300 , 0 . 1367 ) 

- 0.5500 ( 0 . 5500 ) 
( 0 . 1919 , 0 . 2081 ) 

- 0.3762 ( 0 . 3763 ) 
( 0 . 3471 , 0 . 4004 ) 

- 0.1994 ( 0 . 2004 ) 
( 0 . 4996 , 0 . 6015 ) 

- 0.0161 ( 0 . 0598 ) 
( 0 . 6197 , 0 . 8482 ) 

50 

- 0.6833 ( 0 . 6833 ) 
( 0 . 0650 , 0 . 0683 ) 

- 0.6500 ( 0 . 6500 ) 
( 0 . 0960 , 0 . 1040 ) 

- 0.5508 ( 0 . 5508 ) 
( 0 . 1834 , 0 . 2149 ) 

- 0.3889 ( 0 . 3891 ) 
( 0 . 3197 , 0 . 4025 ) 

- 0.0881 ( 0 . 1041 ) 
( 0 . 5414 , 0 . 7823 ) 

100 

- 0.7167 ( 0 . 7167 ) 
( 0 . 0325 , 0 . 0341 ) 

- 0.7000 ( 0 . 7000 ) 
( 0 . 0480 , 0 . 0520 ) 

- 0.6500 ( 0 . 6500 ) 
( 0 . 0921 , 0 . 1080 ) 

- 0.5519 ( 0 . 5519 ) 
( 0 . 1727 , 0 . 2234 ) 

- 0.1983 ( 0 . 2041 ) 
( 0 . 4311 , 0 . 6723 ) 

1000 

- 0.7467 ( 0 . 7467 ) 
( 0 . 0033 , 0 . 0034 ) 

- 0.7450 ( 0 . 7450 ) 
( 0 . 0048 , 0 . 0052 ) 

- 0.7400 ( 0 . 7400 ) 
( 0 . 0093 , 0 . 0107 ) 

- 0.7300 ( 0 . 7300 ) 
( 0 . 0178 , 0 . 0223 ) 

- 0.6500 ( 0 . 6500 ) 
( 0 . 0712 , 0 . 1288 ) 


Table 2: Bias and root mean squared error (RMSE), in brackets, of (u) (first line) and the 95% con¬ 


fidence intervals based on the asymptotic Normal approximation for the blocks estimates (2.13) (second 
line), for a negatively correlated AR(1) process with /3 = 2 and 6 = rj = 0.75. Best values are ticked with 
a plus. 


Example 2.3 (MAR(l) process) Let = ama.x{Xn-i, Cn}) n > 1, 0 < a < 1, denote the max- 
autoregressive process of order one where Xq is a random variable with d.f. Hq, independent of the 
sequence {en}n>i of i.i.d. random variables with unit Frechet distribution. This process, which is a 
special case of the general MARMA(p,q) processes introduced by Davis and Resnick (1989) [1], verifies 
condition D"(un), has extremal index 0 = 1 — a, 0<q:<1 (see Alpuim (1989) [2j and Ferreira (1994) 
m) and upcrossings index r] = 1. 


Tn 

qo.7 

90.8 

90.9 

90.95 

90.99 

5 

- 0.0061 ( 0 . 0092 ) 
( 0 . 9843 , 1 . 0035 ) 

- 0.0039 ( 0 . 0075 ) 
( 0 . 9893 , 1 . 0029 ) 

- 0.0019 ( 0 . 0065 ) 
( 0 . 9947 , 1 . 0016 ) 

- 0.0008 ( 0 . 0056 ) + 
( 0 . 9978 , 1 . 0007 ) 

- 0.0003 ( 0 . 0075 )+ 
( 0 . 9992 , 1 . 0002 ) 

7 

- 0.0143 ( 0 . 0177 ) 
( 0 . 9678 , 1 . 0036 ) 

- 0.0096 ( 0 . 0140 ) 
( 0 . 9759 , 1 . 0048 ) 

- 0.0048 ( 0 . 0110 ) 
( 0 . 9871 , 1 . 0033 ) 

- 0.0024 ( 0 . 0100 ) 
( 0 . 9934 , 1 . 0017 ) 

- 0.0009 ( 0 . 0122 ) 
( 0 . 9982 , 1 . 0001 ) 

10 

- 0.0305 ( 0 . 0341 ) 
( 0 . 9413 , 0 . 9976 ) 

- 0.0202 ( 0 . 0250 ) 
( 0 . 9546 , 1 . 0049 ) 

- 0.0100 ( 0 . 0174 ) 
( 0 . 9736 , 1 . 0064 ) 

- 0.0050 ( 0 . 0148 ) 
( 0 . 9861 , 1 . 0039 ) 

- 0.0013 ( 0 . 0149 ) 
( 0 . 9968 , 1 . 0006 ) 

15 

- 0.0659 ( 0 . 0695 ) 
( 0 . 8935 , 0 . 9747 ) 

- 0.0446 ( 0 . 0494 ) 
( 0 . 9161 , 0 . 9948 ) 

- 0.0232 ( 0 . 0317 ) 
( 0 . 9446 , 1 . 0090 ) 

- 0.0112 ( 0 . 0236 ) 
( 0 . 9702 , 1 . 0074 ) 

- 0.0029 ( 0 . 0214 ) 
( 0 . 9931 , 1 . 0012 ) 

25 

- 0.1495 ( 0 . 1525 ) 
( 0 . 7946 , 0 . 9063 ) 

- 0.1044 ( 0 . 1090 ) 
( 0 . 8378 , 0 . 9534 ) 

- 0.0549 ( 0 . 0636 ) 
( 0 . 8887 , 1 . 0015 ) 

- 0.0280 ( 0 . 0427 ) 
( 0 . 9303 , 1 . 0138 ) 

- 0.0055 ( 0 . 0301 ) 
( 0 . 9866 , 1 . 0025 ) 

50 

- 0.3521 ( 0 . 3539 ) 
( 0 . 5837 , 0 . 7122 ) 

- 0.2624 ( 0 . 2658 ) 
( 0 . 6621 , 0 . 8130 ) 

- 0.1466 ( 0 . 1542 ) 
( 0 . 7653 , 0 . 9416 ) 

- 0.0766 ( 0 . 0924 ) 
( 0 . 8385 , 1 . 0082 ) 

- 0.0153 ( 0 . 0515 ) 
( 0 . 9629 , 1 . 0065 ) 

100 

- 0.6059 ( 0 . 6066 ) 
( 0 . 3447 , 0 . 4435 ) 

- 0.4966 ( 0 . 4982 ) 
( 0 . 4328 , 0 . 5740 ) 

- 0.3123 ( 0 . 3177 ) 
( 0 . 5818 , 0 . 7936 ) 

- 0.1751 ( 0 . 1893 ) 
( 0 . 6988 , 0 . 9510 ) 

- 0.0365 ( 0 . 0838 ) 
( 0 . 9148 , 1 . 0122 ) 

1000 

- 0.9591 ( 0 . 9592 ) 
( 0 . 0363 , 0 . 0454 ) 

- 0.9431 ( 0 . 9431 ) 
( 0 . 0487 , 0 . 0651 ) 

- 0.8932 ( 0 . 8933 ) 
( 0 . 0835 , 0 . 1301 ) 

- 0.7921 ( 0 . 7932 ) 
( 0 . 1420 , 0 . 2738 ) 

- 0.3382 ( 0 . 3843 ) 
( 0 . 4526 , 0 . 8711 ) 


Table 3: Bias and root mean squared error (RMSE), in brackets, of ^n{u) (first line) and the 95% con¬ 


fidence intervals based on the asymptotic Normal approximation for the blocks estimates (2.13) (second 
line), for a MAR(l) process with a = 0.9, 9 = 0.1 and r] = 1. Best values are ticked with a plus. 


Absolute biases tend to be smallest at small block sizes, i.e., Vn = 5, 7, 10 and 15. Furthermore, the 
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results suggest that generally speaking the absolute bias of rj^ (u) increases with but decreases with u 
and the variance of increases with u but decreases with The quality of the estimate depends 

strongly on the choice of the two parameters u and rn, as previously stated. Nevertheless, the blocks 
estimator has a better performance at the 90% and the 95% thresholds. The poor performance at the 
99% threshold is justified by the few observed threshold upcrossings. 

Remark 3 We point out from the previous examples that the stronger (weaker) the dependence between 
upcrossings, i.e. the bigger (smaller) clustering of upcrossings, which implies the smaller (bigger) p, the 
bigger (smaller) values of r^ are required. 

Remark 4 In the estimation of the extremal index the block size rn has commonly been chosen as the 
square root of the sample size, i.e. rn ~ ^/n (Gomes (1993) f77]/ . Ancona-Navarrete and Tawn (2000) 
m, among others). The previous examples show that such a large block size is not a reasonable choice in 
the estimation of the upcrossings index. Choices ranging from Vn = [n/(log n)^'^] to Vn = [n/(log n)^'^] 
seem more reasonable. 


2.1.1.2 Some considerations about the choice of the levels 


In practical applications the deterministic levels Vn previously considered will have to be estimated, 
i.e., they will have to be replaced by random levels suggested by the relation frP[Xi < Vn < X 2 ) ~ n, 
as n —)• + 00 . Nevertheless, these random levels cannot be represented by an appropriate order statistics, 
contrarily to the random levels used in the estimation of the extremal index. 

As in Sebastiao et al. (2013) [23] we shall consider the random levels = A^n-[c„T]:n) used in the 
estimation of the extremal index, more precisely, we shall consider the blocks estimator 

-1 


?n( 


V " IT 
2 ^ 1=1 




Et, 


( +00 ^ 

Y2^j^n{j]Vn) 


(2.14) 


For these random levels Vn, we show in the following result that the blocks estimator r](({vn) is also a 
consistent estimator of p. 


Theorem 2.4 Suppose that for each u > 0 there exists Vn = u^n/c ~ ''^n/c < X 2 ) ~ n 

and -^P{Xi > Vn) ^ t as n ^ +ooj for some r > 0, the conditions of Theorem 2.2 hold for all n > 0, 
then 


VniVn) 


-> rj. 

n^+oo 


Remark 5 If X. has extremal index 0 > 0 and upcrossings index t] > 0, then there exists Vn = u^n/c ~ 
and rju = Or (Ferreira (2006) fUf). 

The high level Vn, considered in Theorem 2.4, must be such that nP(Xi > Vn) = c^r = r„, 
Tn - > +00 and Tn/n - > 0. In the simulation study, presented further ahead, we shall con- 

n^+oo n^+00 

sider the deterministic levels Un = Xn-s-.n, corresponding to the (s -|- l)th top order statistics associated 
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to a random sample {Xi ,..., Xn) of X. The upcrossings index estimator will become in this way a func¬ 
tion o s. In fact, s is replacing = CnT, Cn - > -|-oo. Consistency is attained only if s is intermediate, 

71^+00 

he., s —)■ -|-oo and s = o(n) as n —-|-cx). 

The choice of s, the number of top order statistics to be used in the estimation procedure is a sensitive 
and complex problem in extreme value applications. In Neves et al. (2014) (jl^) an heuristic algorithm 
is used for the choice of the optimal sample fraction in the estimation of the extremal index. The 
optimal sample fraction is chosen by looking at the sample path of the estimator of the extremal index 
as a function of s and finding the largest range of s values associated with stable estimates. Since the 
estimates of the extremal index and the upcrossings index sometimes do not have the largest stability 
region around the target value, this algorithm may give misleading results. 


2.1.2 Short note on other blocks estimators 

If the stationary sequence X also has extremal index 0 > 0 and there exists a sequence of thresholds 


{rxr„}r„>i such that lim r„P(Xi > Ur^) = r and lim r„P(Xi < Ur„ < X 2 ) = then, since (1.4) 

“ rn^+00 r„^+oo 

holds, it follows that 


VrniUrn) — ■ 


log Fr„{Ur„) 


rnP{Xi < Ur„ < X2) r„^+oo 

where Mr„ = max{Xi,..., Xr„} and Fr„{ur„) = P{Mr^ <Ur„). 


> V, 


(2.15) 


The distribution function of the block maximum in (2.15) can be estimated using maxima of 


kn = [n/rn] disjoint blocks or using maxima of re — -|- 1 sliding blocks, as suggested by Robert et al. 

(2009) [2T], in the following ways 


n-Tn + l 


^ '‘'Tl 't- ' n I 

^ 2=1 2=1 




where Mjj = max{Xi,..., Xj}, for 0 < z < j. Robert et al. (2009) [21] showed that Fn {ur„) — 
F^\ur„) has a non-negligible asymptotic variance and is asymptotically uncorrelated with F^{ur„). 
Thus, the sliding blocks estimator is the most efficient convex combination of the disjoint and sliding 
blocks estimators for Fr„{ur„). 

We can now define a disjoint and a sliding blocks estimator for the upcrossings index rj, respectively, 

as ^ 

^j(.. ^ log^n(zZr„) 


and 


C^(WJ:=-- 

VniUrJ ■= - 


log-^n (^r„) 


(2.16) 


(2.17) 


The estimators r]^ and rjfj in (2.16) and (2.17) are not defined when all blocks have at least one 


upcrossing and ^ has already appeared in Sebastiao et al. (2013) 


Remark 6 If for each re > 0 there exists Un'^ and for some r > 0, then P(max{Xi,..., X^} < 


- P{Nn{v>n>) = 0 ) 


(F\ _ 


n^+oo 


> 0 (Ferreira (2006) f^). Therefore, for large re we can replace the 
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probability in the numerator of (2.15) with P{Nr„{ur„) = 0) and estimate it by its empirical counterpart. 
Nevertheless, estimator (2.16) creates synergies in the estimation of the pair {6, rf). We recall that the 
extremal index 9 and the upcrossings index g provide different but complementary information eoncerning 
the oceurrence of rare events. 


Imposing on X a stronger condition than condition A(rt„), involving its maximal correlation coef¬ 
ficients, and considering the threshold sequence {rir„}r„>i deterministic, we have, mutatis mutandis, 
from the arguments used in Robert et al. (2009) [H], that ffH and f](( are consistent estimators of 
g. If in addition, there exists a constant p > 1 such that E Nrf{ur„) = 0(1), as n —+oo, and 


VKigr^UrJ - g) 
that 


-)• 0, then from the central limit theorem for triangular arrays it can be shown 


Tn^+OO 

\ 9niUrJ-g 
where V = [vij] 2 x 2 is a symmetric matrix with 




>+oo 


> A^(0,V) 


uii = - 


V22 = Vi2 = 


1 2 . 

- —-he — 1 

r]u rjv 


- 1 1 n 

2 X —^^— -he — 1 

g^y^ rqy 


and = 




Remark 7 Since V 22 < vn the sliding blocks estimator gf( is more efficient than its disjoint version. 
On the other hand, considering that Cnjkn —>-0, as n ^ + 00 , the blocks estimator g^ is a more efficient 
estimator than g^. 

2.2 Runs declsutering 

Runs declustering assumes that upcrossings belong to the same cluster if they are separated by fewer 
than a certain number of non-uperossings of the threshold. More precisely, if the process X verifies 
condition D^^\un), for some k > 2, of Ferreira (2006) [6], that locally restricts the dependence of the 
sequence but still allows clustering of upcrossings, runs of upcrossings in the same cluster must be 
separated at most by fc — 2 non-uperossings. In Figure 1 we can identify 3 clusters of size 2, since the 
process verifies condition D^^\un). 

Condition D^^\un), k > 2, is said to be satisfied if 

lim nP{Xi <Un < X 2 , N^ k = 0, Nk+i,r„ > 0) = 0, (2.18) 

n^+00 

for some sequence = [n/kn] with {kn}n>i satisfying kn - > + 00 , - > 0, knOni^ ->■ 0, 

n^+oo ^ n^+oo ’ n^+oo 

where an,i„ are the mixing coefficients of the A(tt„) condition, and Nij = Nn{un){[i/n,j/n]) with Nij = 0 
if j < i. 
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This family of local dependence conditions is slightly stronger than D^^\un) of Chernick et al. (1991) 
[3] and (2.18) is implied by 


n ^ PiXi <u< X 2 , Xi<u< Xi+i) 


-> 0 . 

n^+00 


When k 


2 we find the slightly weakened condition D''{un) of Leadbetter and Nandagopalan (1989) 


m 

Under this hierarchy of increasingly weaker mixing conditions D^P{un), k > 2, the upcrossings index 
can be computed as follows: 


Proposition 2.1 (Ferreira (2006) [6j) //X satisfies condition A(un) and for some k >2, condition 
D^^\un'^) holds for some n > 0, then the upcrossings index ofX exists and is equal to rj if and only if 

p(7V3,fc(«(r)) = 0 1 xi < 4") < X2) V 

n^+00 

for each u > 0. 


2.2.1 The runs estimators 

The estimators considered up to now treat each block as one cluster, but under a local dependence 
condition DP'\un)-, k > 3, clusters can be identified in a different way, for example, runs of upcrossings 
or runs of upcrossings separated by at most one non-upcrossings of a certain threshold may define a 
cluster. More precisely, if the process X verifies condition D^^^Un), upcrossing clusters may be simply 
identified asymptotically as runs of consecutive upcrossings and the cluster sizes as run lengths (Ferreira 
(2007) [7j). Therefore, from Propositon 2.1 if for some k > 3 condition D^P{un) holds, the upcrossings 
index can be estimated by the ratio between the total number of A: — 2 non-upcrossings followed by an 
upcrossing and the total number of upcrossings, i.e., the runs estimator given by 

s^n—k _ 

:= ^ ^ 2 . 19 ) 

with Nij = Nn{un){[i/n,j/n]) and = 0 if j < z. 

When k = 3 the runs estimator corresponds to the one fairly studied in Sebastiao et al. (2013) pT| . 
There it was shown that the runs estimator has smaller bias and mean squared error when compared to 
the disjoint blocks estimator rfH. 

^R\k 

Properties such as consistency and asymptotic normality can be proved for the estimators rjn , A: > 4, 
with similar arguments to the ones used in Sebastiao et al. (2013) [23]. Nevertheless, their validation is 
restricted to the validation of a local dependence condition D^^\un), > 4. 

In what follows we propose a way of checking conditions D^P{un), k > 3, preliminary to estimation, 
which turns out to be the only possible solution when dealing with real data. Even when the underlying 
model is known, the verification of conditions DP\un), k > 3, can be cumbersome, so the following 
procedure can in these situations be an important auxiliary tool. A similar approach has been followed 
by Siiveges (2007) [25| to check condition D^'^\un) of Chernick et al. (1991) [3| which is slightly stronger 
than condition D"{un) of Leadbetter and Nandagopalan (1989) [T7] . 
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Condition D^^\un), k > 3, essential for the validity of the runs estimator may be checked by 
calculating 


p^^\u,r) = 


Er=r' h 


{Xi<U<Xi^\, 7Vi + 2,i+i!-l=0, Ni + k,i+r-l>0} 

Y17=i ^{Xi<u<Xi+i_} 


( 2 . 20 ) 


for a sample (Xi,..., Xn) of X, u a high threshold and r the size of the blocks into which the sample is 
partitioned. 

The limit condition D^^\un), for some k > 3, will be satisfied if there exists a path {ui,rj) with 
Ui —)• +00 and rj —?• +oo for which the p^^\ui,rj) —?■ 0. With (2.20) we are looking for the so called 
anti-{un), k > 3, events {iVj+ 2 ,i +2 = 0, > 0 | Xj < u < Xj+i} among the upcrossings for a 

range of thresholds and block sizes. This proportion of anti-D^^\un), k >3, events gives us information 
on how many clusters are misidentified as two or more clusters and therefore information of the bias. If 
such clusters are few, it is plausible to accept the relatively small upward bias despite the possible failure 
of D^^\un), k > 3, and use to estimate r/. This gives us a way to choose, in practical applications, 
among the estimators k > 3, the one to estimate rj. 

In Figures 2 and 3 we plot the proportions p{u,r) = p^^\u,r), considering u the empirical quantiles 
with probabilities ranging from 0.95 to 0.995 and r ranging from 5 to 20, for sequences of size n = 10000, 
of the processes considered in the examples of Section 2, of an AR(2) process and a GARCH(1,1) process 
with Student-t distributed innovations. For the sake of clarity we consider F{u) instead of u in the plots. 

Similar results were obtained for other sample sizes, namely n = 5000 and 30000. For the AR(2) 
process the proportions p^^\u,r) and p^^\u,r) did not differ much from p^^\u,r) which suggest that 
this process also does not verify conditions &^\un) and D^^\un)- The anti-D^^\un)-i k = 3,4, 5,6, 
proportions for the GARCH(1,1) process, plotted in Figure 3, show a more prominent decrease, towards 
zero, for /c = 5 and k = 6. This suggests that conditions D^^\un) and D^‘^\un) are unlikely to hold for 
this process and that k = 5 seems a reasonable choice. 


We recall that the conditions D^^\un), k >3, form a hierarchy of increasingly weaker mixing condi¬ 
tions and therefore, if a process verifies condition D^^'\un) for a fixed k', then it verifies all conditions 
D^^\un) with k > k'. Accordingly we can use any of the estimators , k = k' ,k' -\-l,... ,to estimate ry. 
We simulated for the estimators , A: = 3,4, 5, the mean (E) and the root mean squared error (RMSE) 
from 5000 samples of sizes n = 5000 of the MM process given in Example 2.1. In Figure 4, we find 
the estimated values as functions of s > 1, corresponding to the estimates obtained with [Xn-s-.n)-, 
/c = 3,4, 5, where Xn-s-.n denotes the (s -(- l)th top order statistics. 

DM 

As we can see, from Figure 4, for the MM process the estimator rjn outperforms the estimators rjn 
and rjn , improving the results found in Sebastiao et al. (2013) j2Tj. The estimates have a big stability 
region near the true value ry = 0.5 and the RMSE presents a desirable wide “bathtub” pattern. For the 

D| Q 

processes of Examples 2.2 and 2.3 the best results were obtained with rjn ■ The results show that even if 
we validate condition D^^\un) for some value k > 3 we might obtain better estimates for ry with a runs 
estimator with a higher value of A:. It is therefore convenient to plot sample paths of the runs estimator 
rjn^^[Xn-s:n) for different values of k, as a function of s, and analyze the stability regions. Note however 
that the choice of a too large k may reduce the precision of the new estimators, as can be seen by the 
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Figure 2: Observed proportions p{u,r) = p^^\u,r) for the MM process (top left panel), negatively 
correlated AR(1) process with s = 2 (top right panel), MAR(l) process (bottom left panel) and AR(2) 
process with (pi = 0.93 and (p 2 = —0.86 (bottom right panel). 


values plotted in Figure 4. 


3 Simulation study 


We compare the performance of the several estimators previously presented for a range of different 
processes, having various dependence structures which give both p < 1 and r/ = 1 


In this comparison study we include the estimator that naturally arises from relation (1.4), by repla¬ 


cing 6, T and v with corresponding consistent estimators. Such an estimator has already been considered 
in Sebastiao et al. (2013) |24j . where 0 was estimated with Ferro and Seger’s ([9]) intervals estimator, 
6^^, and r and v by their empirical counterparts. This estimator is formally defined as follows 


= Vn^ (Un) -=0 


:;EI 


where 




1 A 9n{Un) 
1 A 0l{Un) 


'FS 


2=1 




n—1 

2=1 


(3.21) 


max{Ti : l<i<A^ — 1}<2 

max{Tj : l<z<A^ — 1}>2 
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Figure 3: From left to right and top to bottom, observed proportions p{u, r) = r), A: = 3,4, 5 and 

6 for a GARCH(1,1) process with autoregressive parameter a = 0.08, variance parameter /3 = 0.87 and 
innovations 




Figure 4; Mean values (left) and root mean squared errors (right) of the estimators k = 3,4,5, for 
the MM process in Example 2.1 {rj = 0.5). 


with 



2 



N-l 

(iv-1) 


45''-'1 


N-l 


(iv-1) V(r,-i)(r,-2) 


and 0‘^{un) 







denoting Ti, i = 1,. .., N — 1, the observed interexceedance times, i.e., Ti = S’j+i — Si, where 
1 < 5i < ... < Sn < n are the exceedance times. 

The different estimation methods are compared via the estimated bias and variance. Furthermore, to 
analyze the bias-variance trade-off we compare the corresponding root mean squared errors, RMSE. To 
measure the degree of possible misinterpretation caused by considering the variance/standard error as a 
criteria for the quality of the estimator we consider the mean error to standard errors MESE = Bias/\/Var 
(see Frahm et al. (2005) [T0| for further details). 

The values for the statistical quantities considered were obtained with 5000 independent replications 
of the estimation procedures for samples of size n = 5000. Since all estimators of rj depend on a threshold, 
we considered all simulated values for each estimator computed at thresholds corresponding to the 90% 
and 95% empirical quantiles. To reduce sampling variation, for each process all the estimators are 
evaluated using the same simulated data sets. 

For the blocks estimators we compare different block sizes, found to be the best choices from inde¬ 
pendent simulations to those reported here and from the results obtained in Tables 1-3. 

The performance of the estimators for the processes given in Examples 2.1, 2.2 and 2.3 are presented, 
respectively, in Tables 4, 5 and 6. 

Eor a better comparison of the blocks estimators with = 5,10,15 and the runs estimators r/n 
A: = 3,4, which had the best performances in the previous simulations, we plot in Figure 5 the estimated 
mean values and RMSE for levels given by the (s -|- l)th, s = 1, ...,n — 1, top order statistics. 

Some overall conclusions: 

• In general, all estimation methods have small sample bias. For the three processes the blocks 

estimators rj^, = 5, 10,15 and the runs estimators k = 3,4, 5, lead to the smallest absolute 

sample biases. The biggest sample biases are obtained with the sliding blocks estimators and the 
disjoint blocks estimators . The latter has a bit worse performance, except for the AR(1) process 
at the 90% threshold where the sliding blocks estimator and the naive estimator perform worst. 

• The sample variances for all estimators and processes are smaller than 0.05. In the group of the 

smallest observed variance values we find once again the the blocks estimators r/f, = 5, 10,15 

and the runs estimators fyn^^, k = 3,4,5. The conclusions drawn in Remark 7 can be pictured in 
Tables 4, 5 and 6. When ry is on the boundary of the parameter space there is no clear overall 
change in efficiency of the estimators with threshold level. 

• The conclusions previously drawn for the bias and variance carry over to the RMSE, for obvious 

reasons. For the MAR(l) process, with ?y = 1, we obtain the smallest values of the RMSE for all 
estimators, in fact the blocks estimator ?yj^, = 5, and the runs estimator fjn^^ have a RMSE of 

0.006. 

• A large sample bias relative to the sample variance translates into a large MESE, which is true for 

the blocks estimators = 5,10, with the MM and the AR(1) processes. The runs estimators 

present in the majority of the cases values smaller than 1, which indicates that the true upcrossings 
index ry lies within la conhdence band, where cr denotes the standard error. 
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Table 4: Simulation results for the estimators of rj with the MM process of Example 2.1. 


Estimator 

Bias 

Var 

RMSE 

MESE 

Bias 

Var 

RMSE 

MESE 

90.90 

90.95 

Vn (rn 

= 5) 

0.1964- 

0.00033 

0.1972- 

10.8028- 

0.1983- 

0.00069 

0.2000- 

7.5571- 

Vn (rn 

= 10) 

0.0561 

0.00025 

0.0583 

3.5137 

0.0770 

0.00049 

0.0801 

3.4779- 

Vn (rn 

= 15) 

-0.0126+ 

0.00022 

0.0196+ 

0.8459+ 

0.0250 

0.00044 

0.0326 

1.1946 

^ {vn 

= 50) 

0.1195 

0.00404 

0.1353 

1.8797 

0.0717 

0.00228 

0.0862 

1.5002 

^ {Vn 

= 60) 

0.1179 

0.00602 

0.1411 

1.5204 

0.0667 

0.00286- 

0.0855 

1.2471 

^ (r„ 

= 70) 

0.1172 

0.00898- 

0.1507- 

1.2366 

0.0660 

0.00361 

0.0893 

1.0990 

Vn (rn 

= 50) 

0.1156 

0.00202 

0.1241 

2.5719 

0.0705 

0.00129 

0.0791 

1.9663 


= 60) 

0.1124 

0.00283 

0.1243 

2.1139 

0.0660 

0.00164 

0.0774 

1.6305 

Vn (rn 

= 70) 

0.1117 

0.00396 

0.1282 

1.7748 

0.0633 

0.00205 

0.0778 

1.3993 

^R\3 

Vn 


0.0255 

0.00004+ 

0.0262 

4.0038- 

0.0126+ 

0.00004+ 

0.0141+ 

1.9755 

Vn 


-0.0067+ 

0.00009+ 

0.0117+ 

0.6892+ 

-0.0039+ 

0.00009+ 

0.0104+ 

0.4029+ 

^R\5 

Vn 


-0.0376 

0.00013 

0.0394 

3.2787 

-0.0202 

0.00014 

0.0234 

1.7051 

Vn^ 


0.0506 

0.00183 

0.0663 

1.1827 

0.0300 

0.00308- 

0.0631 

0.5415+ 


Best values are ticked with a plus and worst are ticked with a minus. 


• All but the blocks estimator have a better performance the higher the threshold used, since in 
essence they are estimating the non-upcrossing or non-exceedance of the threshold. 

• Estimators obtained from asymptotic characterizations of the upcrossings index have a better per¬ 
formance than estimators obtained from the relation with the extremal index. 

• The good performance of the blocks estimator makes it the best alternative for the runs estima¬ 
tors which need the validation of local dependence conditions. This statement is reinforced by 
Eigure 5, where we see that for an adequate block size the estimated mean of the blocks estimator 
presents a well defined stability region near the true value and the RMSE has a wide “bathtub” 
shape. 
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Table 5: Simulation results for the estimators of rj with the AR(1) process with /3 = 2 of Example 2.2. 


Estimator 

Bias 

Var 

RMSE 

MESE 

Bias 

Var 

RMSE 

MESE 

QO.90 

?0.95 

Vn 

= 5) 

0.0881 

0.00022 

0.0894 

5.9123 

0.1008 

0.00045+ 

0.1030 

4.7764" 

Vn (rn 

= 10) 

-0.0874 

0.00024 

0.0888 

5.6748 

0.0009 

0.00055+ 

0.0235+ 

0.0399 

Vn 

= 15) 

-0.2139 

0.00015+ 

0.2143 

17.6109" 

-0.0776 

0.00055+ 

0.0810 

3.3176" 

^ (rn 

= 50) 

0.1404 

0.00393 

0.1538 

2.2390 

0.1974 

0.01331 

0.2286 

1.7104 

{rn 

= 60) 

-0.0221 

0.00095 

0.0380 

0.7167 

0.2093" 

0.02289" 

0.2582" 

1.3834 

^ (rn 

= 70) 

-0.1428 

0.00020+ 

0.1434 

10.1522 

0.2048 

0.02340" 

0.2556" 

1.3385 

Vn (rn 

= 50) 

0.3980" 

0.04007" 

0.4455" 

1.9883 

0.1890 

0.00698 

0.2066 

2.2622 

Vn {fn 

= 60) 

0.3100" 

0.02936" 

0.3542" 

1.8094 

0.1949 

0.01014 

0.2194 

1.9349 

Vn 

= 70) 

0.1953 

0.01983 

0.2408 

1.3870 

0.2050" 

0.01545 

0.2397 

1.6490 

^R\3 

Vn 


0.0006+ 

0.00039 

0.0198+ 

0.0279+ 

0.0007 

0.00077 

0.0278+ 

0.0255 

^14 

Vn 


0.0004+ 

0.00039 

0.0198+ 

0.0193+ 

0.0006+ 

0.00077 

0.0278+ 

0.0208+ 

Vn 


-0.0622 

0.00044 

0.0657 

2.9574 

0.0004+ 

0.00077 

0.0278+ 

0.0159+ 

Vn^ 


0.2424 

0.00040 

0.2432 

12.0808" 

0.1566 

0.00394 

0.1687 

2.4965 


Best values are ticked with a plus and worst are ticked with a minus. 


4 Application to financial data 

We consider the performance of the estimators under study which have previously better performed, 
i.e. the blocks estimator rj^ and the runs estimators k > 3, when applied to the analysis of the 

log-returns associated with the daily closing prices of the German stock market index DAX, collected 
from 3 January 2000 to 29 June 2012. This series was analysed previously by Sebastiao et al. (2013) [M] 
and is plotted in Figure 6: DAX daily closing prices over the mentioned period, xt, and the log-returns, 
100 X (Inxi — lnxi_i), the data to be analyzed, after removing null log-returns (n = 3177). 

As stated in Sebastiao et al. (2013) [23|, Klar et al. (2012) [20j have analyzed the DAX German 
stock market index time series and concluded that the GARCH (1,1) process with Student-t distributed 
innovations is a good model to describe these data. The simulated values of the anti-D*^^^ (un), k = 3,4,5, 6 
proportions found in Figure 3 for the GARCH(1,1) process and the conclusions there taken led us to 
consider the runs estimators with = 4, 5 and 6. For the blocks estimator, the choice of r„ = 5, 7 
and 10 seemed adequate for the sample size considered. 
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Table 6: Simulation results for the estimators of r/ with the MAR(l) process with a = 0.9 of Example 

2.3. 


Estimator 

Bias 

Var 

RMSE 

MESE 

Bias 

Var 

RMSE 

MESE 

90.90 

90.95 

Vn (rn 

= 5) 

-0.0019+ 

0.00004+ 

0.0065+ 

0.3006+ 

-0.0008+ 

0.00003+ 

0.0056+ 

0.1420+ 

Vn (rn 

= 10) 

-0.0100 

0.00020 

0.0174 

0.7001 

-0.0050 

0.00019 

0.0148 

0.3583 

Vn (rn 

= 15) 

-0.0232 

0.00046 

0.0317 

1.0804 

-0.0112 

0.00043 

0.0236 

0.5397 

^ {Vn 

= 50) 

0.3304- 

0.01156 

0.3474- 

3.0723 

0.2586 

0.01341 

0.2834 

2.2335 

^ {Vn 

= 60) 

0.2986 

0.01363 

0.3206 

2.5580 

0.2248 

0.01387 

0.2538- 

1.9086 

^ (r„ 

= 70) 

0.2759 

0.01526- 

0.3023 

2.2332 

0.2066 

0.01429- 

0.2386 

1.7283 

Vn (rn 

= 50) 

0.3270- 

0.00591 

0.3359- 

4.2524- 

0.2552- 

0.00594 

0.2666- 

3.3119- 

Vn (rn 

= 60) 

0.2949 

0.00687 

0.3063 

3.5595- 

0.2241 

0.00626 

0.2377 

2.8332 

Vn (rn 

= 70) 

0.2726 

0.00806 

0.2870 

3.0372 

0.2028 

0.00695 

0.2192 

2.4321 

^R\3 

Vn 


-0.0013+ 

0.00003+ 

0.0055+ 

0.2518+ 

-0.0009+ 

0.00003+ 

0.0060+ 

0.1451+ 

^R\4 

Vn 


-0.0035 

0.00007 

0.0091 

0.4140 

-0.0021 

0.00009+ 

0.0095 

0.2250 

^R\5 

Vn 


-0.0066 

0.00014 

0.0135 

0.5666 

-0.0036 

0.00014 

0.0125 

0.3019 

Vn^ 


0.0324 

0.01645- 

0.1323 

0.2523 

0.0426 

0.02779- 

0.1720 

0.2558 


Best values are ticked with a plus and worst are ticked with a minus. 


Clusters of upcrossings of high returns can bring serious risk to financial investors so estimates of 
the upcrossings index for the log-returns are plotted against s > 1 in Figure 7, in a linear scale and in 
logarithmic scale to better observe estimates at high levels. Note that Figure 7 also enhances the effect 
that the block size has on the estimates. 

The sample paths of all estimators have a stability region around the value rj = 0.8, which can be 
identified in the plot of the estimates in a logarithmic scale of Figure 7. The longest stability region is 
obtained with the blocks estimator for = 7. The higher value of rj obtained with the runs estimator 
comparatively to the ones obtained with the runs estimators and can be justified by the 
fact that condition does not hold and therefore runs of upcrossings are over evaluated. 

A value of t] signihcatively less than 1 shows certain short range dependence reflected in the clustering 
of high level upcrossings. This may be interpreted as a pattern in the occurrence of the extreme values, 
that is, once a large profit in the asset return has occurred it will be followed by a loss (negative return) 
and we can expect a period of large prohts followed by losses (values crossing some threshold). The 
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Figure 5: From left to right, mean values (top) and root mean squared errors (bottom) of the estimators 
with Vn = 5,10,15 and A: = 3,4, for the MM process in Example 2.1 (77 = 0.5), the AR(1) process 
with /3 = 2 in Example 2.2 (77 = 0.75) and the MAR(l) process with a = 0.9 in Example 2.3 (77 = 1). 
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Figure 6 : DAX daily closing prices (left) and daily log-returns (right), of DAX from 2000 to 2012, with 
3177 observations (successive equal prices excluded). 

average length of this period being the inverse of the upcrossings index. 
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Figure 7: Estimates of the upcrossings index for the daily DAX log returns plotted against s > 1, in a 
linear scale (left) and in a logarithmic scale (right). 


5 Conclusions 

New blocks estimators for the upcrossings index tj, a measure of the clustering of upcrossings of high 
thresholds, have been proposed in this paper. These estimators have the advantage of not requiring the 
validation of a local dependence condition, contrarily to the known runs estimator of the upcrossings 
index. The study of the properties of these blocks estimators, namely their consistency and asymptotic 
normality, has been carried out in Section 2. The family of runs estimators of t] has been enlarged and an 
empirical way of checking local dependence conditions that control the clustering of upcrossings provided. 

The estimators proposed in this paper improve substantially on existing methods of estimating the 
upcrossings index. Nevertheless, the simulations performed, in Sections 3 and 4, tell us that there is still 
space for improvements. In future, a discussion on the optimal choice of the key parameter has to be 
considered as well as bias reduction techniques, namely the generalized Jackknife methodology. 


Appendix A: Proofs for Section 2 


A.l Proof of Theorem 12.II 

Let us consider the sample Xi ,..., divided into knjcn disjoint blocks. We can then apply the 

arguments used in Lemma 2.1 of Ferreira (2006) [6] and conclude that 


P(iV[n/c„](^n) = 0) ~ = 0)- 

Now, from the definition of the upcrossings index rj we have P(A[„/c„] ('^n) = 0) 
applying logarithms on both sides of (A.l), (2.9) follows immediately. 


n^+oo 


(A.l) 
> e~^'^, hence 


For the conditional mean number of upcrossings in each block, we have from (2.9) and the definition 
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of the thresholds Vn 


I Nr„{Vn) >0] — 


^ — jl ^Tni'^n) > 0) 

Pi^rniyn) > 0 ) 

E[Nr^{Vn)\ _ rnP[Xi < Vn < X2) 


PiNr„{Vn) > 0) P{Nr„{Vn) > 0) 

TnCny/n 

CnTJl^/kn’ 


which completes the proof since ~ n. 


□ 


A.2 Proof of Theorem 12.21 

It suffices to show that 


+ 00 


lim y^jTTn{j]Vn) 

-A-hno ^ ^ 


n^+oo 


= ^ 




i=i 


and 


+ 00 


j{'^n{j]Vn) -7^n{j]Vn)) - 0 . 

n^H-CO 


(A.2) follows immediately from Theorem 2.1 since 

+ 0O 

lim y^innij;Vn) = hm E[Nr^{Vn) \ Xr„{Vn) > 0] = r]~^. 

n —>■+00 ^ ^ n—>-l-on 


i=i 


n^+oo 


(A.2) 


(A.3) 


To show (A.3), lets start by noting that Theorem 2.2 of Hsing (1991) [T3] holds for nn{j',Vn) and 
^n{j',Vn), if Yni, the number of exceedances of Vn in a block of size r^, is replaced by Nr„{vn)- The proof 
now follows from considering a = 1 and T{j) = in this theorem for Nr„{vn) and verifying that 

conditions (a), (c) and (d) of this theorem follow from the assumptions of Theorem 2.1, (2) and (3), (b) 
follows from Theorem 2.1 and (e) from the fact that 

hm —E[Nr„{vn)] = lim —rnP{X<Vn < X 2 ) = v. 

n—^+00 Cn n^+00 Cn 

□ 


A.3 Proof of Theorem 12.31 


From Cramer-Wald’s Theorem we know that a necessary and sufficient condition for (2.10) to hold is 
that for any a, 6 £ H 

/ \ 

^-n'^ ( «E (E®K) - E[N^}{Vn)]) + 


2 = 1 


2=1 




n ^+00 


> N{Q, ^{a^rja^+ 2ab+ b^r])). 


(A.4) 


We shall therefore prove (A.4). For this, lets consider 


'i'P n 

^rl-lS'^ri) = E 

j=(i-l)rn+l 




(0 


23 
















and note that 


2=1 




a (iV«K) - E[N(^{vn)]) + & E 

2=1 

hfi kfi 


2 = 1 
hn 


2=1 


2=1 

kn 




rn — l 


{K:.^jv.)>oy 


2=1 


rr(*) 0,_ ^—nl '^[^/rr(®) 


K)>0, K)=o} ^ K)>o. K)=0}J J • 


(A.5) 


Now, since (A.5) holds, to show (A.4) we have to prove the following 

/ kfi kn 


.-1/2 


E {"S-lM - EKl,JVn)\) + 4E “ 


2 = 1 


2=1 








n^+oo 


in ^ ^ In ^ r).-^-kno 


> AA(0, v{a^r](T‘^+ 2ab+ b^ri)), 


(i). 


2 = 1 

kn. 


.-1/2 


Eh 


2 = 1 


n^+oo 


NAfS>0, <Li„(«n)=0} f/«_,„(t;„)=0}jy n^+oo 


> 0 . 


Lets first prove (A.6). The summands in Yli=ii'^Ur^-i i'^n) + ^2, (i) 




(A.6) 

(A.7) 

(A.8) 

are functions of 


indicator variables that are at least In time units apart from each other, therefore for each f G IR 

kn 


E 


exp 


f^n / 

^ (aUl!'ES^n) + M 




(A.9) 


2=1 


Hl; \tc-^/‘^[aU^l,^{vn) + b^ 


{U^^.lAvn)>0} 


E 16fcn,crn,z„) 


(A.IO) 


where, i is the imaginary unit, from repeatedly using a result in Volkonski and Rozanov (1959) [26] and 


the triangular inequality. Now, since condition A{vn) holds for X (A.IO) tends to zero and so we can 
assume that the summands are i.i.d.. Therefore, in order to apply Lindberg’s Central Limit Theorem we 
need to verify that 

2" 


kn ^ 
—E 

Cn. 


®i7r„-Z„(Un) + 


since ^E'^[Ur„-iSvn)] 


n^+oo 


n^+oo 

> 0, and Lindberg’s Condition 


> u{a^ria‘^ + 2ab + rrj), 


(A.ll) 


kn ^ 
—E 


(^aUr„-l„iVn) + ^{(af/,.„_i„(i;n)+M{c/^^_,^(t,„)>o})^>cc„} 


n^+oo 


> 0, for all e > 0, 
(A.12) 
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with Ur„-l^{Vn) = 


From the definition of Nrl^{vn) and and assumption (2) we have that 

hn 


/ 




(A. 13) 


n^+oo 


> 0 , 


with Vr„{vn) = Now, by Cauchy-Schwarz’s inequality 

^E[Nr„{Vn)VrM] < E[Nl{Vn)]E[Vl{Vn)] 

thus 

^E[Ul_i^{Vn)] = ^E[{NrM - VrMf] 

Cn Cn n^+oo 

On the other hand, since kn{rn — In) ~ n, Theorem 2.1 implies that 


kn ^ 
—E 

Cn. 


I 


{Vrn-ln (^n)>0} 


= —P{Ur^-l^{Vn) > 0) -^ r]V. 

Cri n^+oo 


(A.14) 


(A.15) 


Furthermore, by definition (2.8) 
k 


—E 


-I Cn Cn 


n^+oo 


jy. 


(A. 16 ) 

(A.14)-(A.I 6 ) prove (A.11) and since Lindberg’s Condition follows immediately from assumption (1), 


(A.6) is proven. 


Finally, since Theorem 2.1 of Hsing (1991) [l5] holds for Nr^{vn), it implies (A.7) and (A.8) because 
^E[Vi^Vn)] —^ 0 by (|aT^ and 


n^+cxo 


—E[n^ 


k 


{N^'^iv„)>o, u^;>ijv„)=oy c, 


rii) 


-nr] = '^{PiN^H^n) > 0) - P(t/« (u„)) > 0) 0 


n^+oo 


by Theorem 2.1 and (A.15). This concludes the proof. 


□ 


A.4 Proof of Corollary 2.1 


Since -)■ rj, by Theorem 2.1, ^E[Nr„{vn)] — 

n—¥-\-oo n^+oo 


> u, and c.^^Yli=iiE^rl!ivn)-E[Nyj{vn)]) — 

n —>-+00 




r(*)/ 


0, by Theorem 2.1 of Hsing (1991) [15] which holds for NrJ'P{vn), the result now follows from the fact 
that 

1 f 


V^(j]n -^n) = — 




2 = 1 


\0r)ivn)>O} ■^[^{7V«(r;„)>0} 


-VnCyP^ Y. {N^i^iVn) - (^^n)]) 


2=1 


and Theorem 2.3. 
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A.5 Proof of Theorem 12.41 


Since (A.2) holds, we only need to show that 

H-OO 


{T^n{j]Vn) - vrn(j>n)) -^1 0. 

n^+oo 


(A.17) 


Lets start by noting that for Vn'^ such that P{Xi > Vn'^) ~ CnT/n as n ^ +oo and e > 0 we have 
kv 


lim lim — 

e^O n—^+oo Cn 


E 






-E 






= lim lim 


kriTri 


€—^0 Ti —^H-CXD C72 

|. ^nXr 

= lim lim 


e^On ^+00 Cn 


(P(Ai < < A2) - P(Xi < 4 ") < A2)) 

(P(Ai > r;("+^)) - P(Xi > 4"))) = 0. 


(A.18) 


(A.18) proves condition b) of Theorem 2.3 in Hsing (1991) [IS] which holds for Nj.„{vn), where T{j) = 
jl[|j>i}. The other conditions have been verified in the proof of Theorem 2.2 as well as the conditions of 
Corollary 2.4 in Hsing (1991) [T5| for Nr„(vn)- Therefore (A.17) holds, completing the proof. □ 
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